Phase switching in a voltage-biased Aharonov-Bohm interferometer 
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Recent experiment [Sigrist et al., Phys. Rev. Lett. 98, 036805 (2007)] reported switches between 
and TT in the phase of Aharonov-Bohm oscillations of the two-terminal differential conductance 
through a two-dot ring with increasing voltage bias. Using a simple model, where one of the dots 
contains multiple interacting levels, these findings are explained as a result of transport through 
the interferometer being dominated at different biases by quantum dot levels of different "parity" 
(i.e. the sign of the overlap integral between the dot state and the states in the leads). The 
redistribution of electron population between different levels with bias leads to the fact that the 
number of switching events is not necessarily equal to the number of dot levels, in agreement with 
experiment. For the same reason switching does not always imply that the parity of levels is strictly 
alternating. Lastly, it is demonstrated that the correlation between the first switching of the phase 
and the onset of the inelastic cotunneling, as well as the sharp (rather than gradual) change of 
phase when switching occurs, give reason to think that the present interpretation of the experiment 
is preferable to the one based on electrostatic AB effect. 

PACS numbers: 73.23.-b, 73.23.Hk, 73.63.Kv 



I. INTRODUCTION 

The phase of a transmission coefficient through a 
quantum dot (QD) can be experimentally studied us- 
ing an Aharonov-Bohm (AB) interferometer with the 
QD embedded in one of its arms,^'^ Since the early 
experiments^, large body of theoretical work has been 
devoted to study of the rich phase behavior experimen- 
tally discovered. When studied in two terminal geom- 
etry such an interferometer exhibits " phase rigidity" , 
i.e. the AB oscillations of conductance have either 
maximum or minimum at zero magnetic field, which 
corresponds to a transmission phase equal to or tt. 
Phase rigidity has its origins in time-reversal symmetry 
(i.e. unitarity of the transmission matrix).- In "open" 
interferometers^ or two-terminal interferometers in non 
equilibrium conditions^ the phase rigidity may be bro- 
ken, i.e. the phase may change continuously. Particularly 
interesting are phenomena connected to " phase lapses" , 
which are abrupt changes of the transmission phase as 
a function of the plunger voltage in the Coulomb block- 
ade valley (i.e. between two successive Coulomb blockade 
resonances). While there has not been yet general agree- 
ment about the physics underlying this phenomenon^^i^ 
it is believed to be the consequence of the system un- 
dergoing transitions between regimes where transport is 
dominated by different levels 

Recently, an interesting new phase phenomenon has 
been observed"'^^ in an AB interferometer under non- 
equilibrium conditions (i.e. with finite bias applied be- 
tween the source and the drain). In this experiment the 
differential conductance of a two terminal interferome- 
ter was measured as a function of the bias voltage and 
magnetic flux through the interferometer. Both arms of 
the interferometer contained quantum dots, which were 



tuned to the Coulomb blockade regime, thus allowing 
only for cotunneling transport. As a function of bias sev- 
eral switches of the phase of the AB oscillations between 
the values and tt were observed with the first switch 
coinciding with the onset of inelastic cotunneling in one 
of the dots. 

A possible explanation for this effect may be as the 
result of electrostatic AB oscillations, as was studied in 
Refsjlliil^,. Some features of the experiment, however, 
are inconsistent with such an explanation. In particu- 
lar, as mentioned above, there is apparent correlation 
between the onset of the phase switching and the onset 
of inelastic cotunneling. In addition, the phase of AB 
oscillation, which at finite bias is not limited by symme- 
try to or vr, is in fact changing abruptly between these 
values as a function of bias, whereas in Ref 12 it changed 
smoothly. 

An alternative explanation (motivated by " population 
switching" studiesS'^) is that, depending on the bias volt- 
age, the elastic contribution to the transport through 
the quantum dot may be dominated by different orbital 
states within the dot. Depending on the sign of the ma- 
trix elements between the state in the dot and the states 
in the leads ("parity" of the state), an electron can ac- 
quire phase or tt, additional to the other phases ac- 
quired when traversing the interferometer. This results 
in AB oscillations having maximum or minimum depend- 
ing on the parity of the current-transmitting state. Thus, 
if the two adjacent QD states have different parity, the 
phase switching may occur in the regimes when the state 
dominating the coherent transport changes from one to 
the other, which may happen with increasing bias. 

From the simple theory described above one would 
naively expect that, if more than two levels in the dot are 
energetically accessible, multiple switching events would 
occur only when the levels have alternating parity, and 
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the number of switching events would then be equal to 
the number of parity changes. This presumption contra- 
dicts the experiment where the number of switchings can, 
in fact, be greater than the number of inelastic onsets, 
i.e. greater than the number of the QD levels involved. 
In what follows we will show that neither of these expec- 
tations is quite true, and a simple model can reproduce 
all the features of the experiments^ as result of redistri- 
bution of electron population between QD states. We 
also point out that at finite bias "population inversion" 
is not a necessary condition for the phase switching to be 
observed. 

The plan of the paper is as follows: in section ITTl we 
first introduce our model of the interferometer, which, 
after a Schrieffer- Wolff transformation,-- allows to study 
cotunneling processes in terms of a single particle scat- 
tered from a local pseudo-spin degree of freedom. In 
section Hm this model is solved in terms of rate equations 
and the differential conductance is obtained as a function 
of the Aharonov-Bohm flux threading the interferome- 
ter and source-drain bias. We compare results with the 
experiment, showing, in particular, that the number of 
switching events can be greater or less then the number 
of levels, whose energy spacing is within the bias window. 
We discuss how the switching behavior may occur even 
when the parity of levels is not strictly alternating, and 
compare our results with the electrostatic AB scenario. 

II. MODEL 



of computational convenience. 

In the present experimental context, the Kondo effect 
is irrelevant and its effects will be neglected in the present 
approach. In addition, for computational simplicity, we 
assume that each level carries a different quantum num- 
ber, which are not mixed in the leads (i.e. we disregard 
difference between spin and orbital channels). The re- 
sults obtained below can be shown to be correct also in 
the case when the number of channels in the leads is not 
equal to the number of states in the quantum dot (e.g. 
when all the levels in the dot are coupled to the same lead 
channel). However, since the dot levels in such a system 
are connected via tunneling into the leads, the correct 
calculation is more complicated and requires a decompo- 
sition procedure similar to that developed in Refill. 




FIG. 1: (color online) Schematic diagram of the modeled de- 
vice. Two infinite tight-binding chains, possibly with different 
chemical potential are coupled to an Aharonov-Bohm inter- 
ferometer, one of the arms of which contains a multi-level 
quantum dot. 



A. Hamiltonian 

We describe the interferometer by a tight-binding 
model, schematically shown in Fig. [TJ The Coulomb 
interaction is present in one of the dots, which contains 
several electron orbitals (interacting arm), whereas the 
dot in other arm of the interferometer (reference arm) 
has only one level. The chemical potential is set out of 
resonance, so that the transport via the interacting arm 
occur only by means of cotunneling processes. The refer- 
ence arm can be characterized by its one-particle trans- 
mission coefficient, and placing a level in it is a matter 



Cotunneling processes can be described as second order 
transitions between different states which contain exactly 
one electron one in the doti^ By performing a Schrieffer- 
Wolff transformatior^i^ (Appendix |^ we can reduce the 
problem to that of one particle scattered by local pseudo- 
spin degree of freedom. As a result, we can describe the 
system by the following (one-particle) Hamiltonian: 

n^HL+HR + HD+Wref+W2 + Wi + WQ, (1) 

where 
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with ric = d^da). Here a,/? refer to the different dot 
levels/lead channels, whereas n is the tight-binding site 
label with n > 0(n < 0) in the right (left) lead. c+ „(ca,„) 
creates (annihilates) an electron in channel a on site n, 
whereas d^{da) creates (annihilates) electron in state a 
in the quantum dot with energy . Vl , Vr are the po- 
tentials applied to the leads, related to the chemical po- 
tentials as hl,r = cf + Vl.r (ep is the Fermi energy), t 
is the hopping integral. 

The terms H^, describe non-interacting leads, Hd 
is the Hamiltonian of the dot, W2 and W% describe elastic 
cotunneling through the quantum dot when it is empty 
or occupied by one electron respectively. W4 describes 
inelastic cotunneling, i.e. the processes when an electron 
incident in channel a is scattered into channel (3, while 
the dot changes its state from /3 to a. Wref describes tun- 
neling through the reference arm, which we have taken 
as channel- independent. The AB phase, 4>ab, can be 
attached to the matrix elements for the reference arm: 
LOri = uj*j. = \ujri \ exp {i(j)AB)- The explicit dependence of 
the tunneling parameters v and oj on the original param- 
eters appear in Appendix VK\ 



where |0) is the state with no particles. The wave vectors 
qp , kf3 are defined as solutions of the equations 

E ^ + Vl - 2tcosqf3, 

E = + VR-2t cos kp (5) 
(el (9) Vl + f-q and eR{k) = Vr, + eu, where ej, = 
— 2icosfc are electron tight-binding eigenenergies in the 
two leads). Let M be the number of the channels, then 
Ai,A/?,Bi,Bfl are vectors with components and 
the amplitude scattering matrix can be defined as (more 
details are given in Appendix IB|) 

In order to normalize the wave function and the scatter- 
ing matrix to unit incident flux, we multiply S by velocity 
matrices: 



B. Scattering matrix 



where 



The energy of the system "QD+incident electron" is 
conserved, so we can start with calculating the scattering 
matrix, which we do by solving the Schroedinger equa- 
tion. 



Sa.a'Sp,0''2t sing^, 
: Sa^a'Si3^f3'2t sin kp. 



(8) 



H\^)^E\-^), 



(3) 



Transition rates 



where E is the energy of the system "dot-|-incident elec- 
tron" for the wave function written as 



a, 13 Ln= — 00 
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a/3 -ife^r 
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In order to get the rate equations we write down the 
rate at which electrons incident from channel a of lead 
fi{— L, R) are being scattered into channel a' of lead fj,', 
whereas the impurity changes its state from (3 to (3' . (In 
fact, only two types of processes are possible: a) elastic 
scattering, when a' = a,P' = (3, and b) inelastic scatter- 
ing, when a' = (3,(3' = a {a ^ (3).) 
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ep,). (9) 
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Here we wrote explicitly the dependence of the scatter- 
ing matrix on the energy of the system " QD-t- incident 
electron", whereas /o(e) = + exp {{e — €F)/{kBT))] 
is the Fermi distribution function; Vl — Vj^ = V is the 
bias voltage. 

Level occupation numbers. The resulting rate equa- 
tions, describing transitions between different states of 
the quantum dot are: 



(10) 



w 



,/3;/3 



' = EEw^^^"''' = E<J''''/3^/3', (11) 



where Pfiit) is the probability that the electron occupy- 
ing the quantum dot is in the state with energy ep. In 
a stationary state the occupation numbers, Pp can be 
calculated from equations 



dt 



Pp{t) = 



(12) 



supplemented by the normalization condition for the to- 
tal number of electrons occupying the dot. 



(13) 



Current. The lead-to-lead current is given by / = 

^coh H" lincoh •, 



Q,/3 



^incoh 



(14b) 



a, 13 



where we have explicitly separated the elastic and inelas- 
tic contributions. 

The differential conductance is defined as 



G{V,(t> 



AB 



dljV, (t>AB) 

dV 



(15) 



where (j)AB — 27r$ / $o is the Aharonov-Bohm phase (pro- 
portional to the magnetic flux, $, divided by the mag- 
netic flux quantum, $0 = hc/e). 



III. RESULTS 

We provide results obtained from the above equations 
for an interferometer that has three levels in the inter- 
acting quantum dot. For the sake of clarity we chose 
£3 — ^2 = 2(e2 — ei), so that inelastic cotunncling pro- 
cesses involving levels 1 and 2 begin at lower bias than 
the cotunneling processes that involve levels 2 and 3. The 
couplings of the levels to the leads increase with their 
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FIG. 2: (color online) Normalized oscillating part of the dif- 
ferential conductance, Eq. (|16|l , as a function of magnetic flux 
and bias voltage for the set of parameters appearing in Ap- 
pendix [C] The conductance maximum as a function of flux, 
as zero bias, becomes a minimum and then a maximum again, 
as bias is increased. 



energies, r2/ri ~ 3.2, Fa/Fi ~ 9. The complete set of 
parameters appear in Appendix [C] 

Our main result is Fig. [21 which shows a plot of the 
oscillating component of the differential conductance as a 
function of the magnetic flux through the AB ring (hor- 
izontal axis) and the bias voltage applied to the inter- 
ferometer (vertical axis). For comparison with the ex- 
perimental results we normalized the conductance in the 
same way as in Reffiol , i. e. the quantity shown is 



g{V,(t)AB) = 



G{V,^ab)-G 

mean iv) 



(16) 



where Gmean (V) is the differential conductance averaged 
over the period of the magnetic flux, 4>ab- Comparing 
this figure to Fig. 3 of the experiment by Sigrist et al.fiS 
we observe the same characteristic features: (i) depend- 
ing on the bias voltage, AB oscillations of the differential 
conductance may have maximum or minimum at zero 
magnetic flux; (ii) the amplitude of the conductance os- 
cillations relative to the non-oscillating background de- 
creases with increasing bias; and (iii) the number of 
switching events (i.e. instances when conductance AB 
oscillations change from maximum at 4>ab = to mini- 
mum or vice versa) is not equal to the number of inelastic 
onsets (in the case at hand we have four switchings, but 
only three inelastic onsets, corresponding to e2 — ci, £3 — £2 
and €3 — El, respectively). Let us discuss these features 
in greater detail. 

The mean value and the oscillating component (at zero 
magnetic flux) of the differential conductance as func- 
tions of bias are shown in Figures [3] and [3] respectively. 
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FIG. 3: (color online) Differential conductance (solid blue), 
its derivative (dash green) and AB phase of conductance os- 
cillations as functions of bias voltage (dotted black). The 
onset of the inelastic processes, corresponding to V = ej — 
(arrows) are evident. 



FIG. 4: (color online) The oscillating contribution to the dif- 
ferential conductance at zero field (solid blue) and AB phase 
of conductance oscillations (dash green) as a function of bias 
voltage. 



In Fig. [3] one can see clearly the three inelastic onsets, 
that is three sharp increases of conductance at bias volt- 
age equal to spacing between any two of the three levels 
in the system (Arrows corresponding toV = ti — tj). For 
clarity we show on the same plot the first derivative of the 
conductance and the phase of conductance oscillations 
(both multiplied by constant factors). The phase is de- 
fined as the sign of the AB oscillations relative to Gmean, 
i.e. it is TT if the oscillations have maximum around zero 
magnetic flux and otherwise. It is necessary to point 
out that at finite bias the Onsager-Biittiker symmetries 
for a two-terminal structure do not require that the phase 
should be restricted to 0, tt,^ therefore our definition of 
phase is not mathematically strict. Nevertheless, for the 
set of parameters used here (see Fig. [5]) an extremum is 
found at zero phase even at finite voltage. 

The oscillating component of the conductance (Fig. [4|) 
changes its sign four times. The first change occurs when 
the bias reaches value V = £2 — ei, i-e. when inelastic 
cotunneling process becomes possible, in which case an 
electron occupying level 1 leaves the dot, whereas an- 
other electron can enter the dot and occupies level 2. 
While such processes do not contribute to the coherent 
AB oscillations (they change the state of the dot, thus 
leaving "trace in the environment"—), they play an im- 
portant role in changing the occupation of the different 
levels of the dot. In particular, population of level 2, 
which was for small bias occupied only due to finite tem- 
perature, increases, as can be seen in Fig. [SJ where we 
show how the populations of all three levels in the inter- 
acting arm of the interferometer depends on the applied 
bias (at (/)a_b = 0). AB oscillations of conductance re- 
sult from elastic cotunneling processes, in which electron 
can tunnel through any of the dot levels. The relative 
weights of these processes are proportional to the proba- 
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FIG. 5: (color online) Populations of QD levels as functions 
of bias voltage. With increasing voltage, and with the onset 
of inelastic processes (arrows) there is a significant change in 
the relative occupations of the different levels. 

bility for the respective level to be occupied, and there- 
fore for low bias all tunneling processes are suppressed, 
except those that happen via level 1 or via intermediate 
state in which the dot is occupied by two electrons (the 
latter will be suppressed, if the Coulomb interaction is 
large). Thus, the increase of the occupation of level 2 
enhances the weight of the corresponding elastic cotun- 
neling processes, whose contribution to AB oscillations 
has phase opposite to that of level 1 (due to the opposite 
parity of these two levels). Due to stronger coupling of 
level 2 to the leads the oscillations of opposite phase even- 
tually outweigh those due to level 1 and phase switching 
occurs. 

This interplay between the level coupling strength and 
its population leads to some interesting effects that hap- 
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pen when the bias is further increased. At the second 
inelastic onset, i.e. when V — €3 — 62, occupation of level 
3 starts to grow due to transitions from level 2. However, 
the bias is still not high enough to excite electrons from 
level 1 to level 3 directly, therefore level 3 is still being 
intensively depopulated by transitions to level 1, which 
results in simultaneous increasing of the population of 
levels 1 and 3, and saturation (or even a decrease) of the 
population of level 2. Therefore the second switching of 
AB phase occurs due to the change in the relative popu- 
lations of levels 1 and 2. It is necessary to stress the role 
of the coupling strengths on the population redistribu- 
tion: for low temperatures the population of level 3 may 
be extremely small, but the switching still occurs since 
level 2 is effectively depopulated by transitions to level 1 
via level 3. 

Further increase of the bias leads to increase of the 
population of level 2, which again outweighs that of levels 
1 and 3 and causes the third switching. Finally, at the 
biases greater than 63 — 61, direct population of level 3 
from the lowest energy level 1 starts and the contribution 
of transport with level 3 occupied eventually takes over 
all other contributions and results in the fourth switching 
occurrence. 

The contributions of different levels to the AB oscilla- 
tions are shown in Fig. ]6\ Let us also point out that a 
contribution of a level to AB oscillations may not always 
have the same sign. Thus, there is a regime in which two 
levels of different parity (levels 2 and 3) give the same sign 
contribution to AB oscillations. The reason for that is 
finite value of Coulomb interaction, which allows for two 
cotunneling processes: a) hole tunneling via the occupied 
level (with intermediate state being an empty dot), and 
b) electron tunneling via an unoccupied level (with inter- 
mediate state being the dot twice occupied). These two 
processes have energy denominators of different signs. In 
other words, the former process involve exchanging the 
two electrons and therefore contributes with the sign dif- 
ferent from the latter process. Thus, for example, the 
contribution of level 3 to the AB oscillations is negative 
as long as this level is unoccupied, but it becomes posi- 
tive when the population of this level becomes non-zero. 
This change of sign does not occur in the limit of large 
Coulomb interaction, when only hole tunneling (or only 
electron tunneling) is possible. 

Coherent (elastic cotunneling) and incoherent (inelas- 
tic cotunneling) contributions to the conductance are de- 
picted in Fig. [7l Inelastic cotunneling never contributes 
to AB oscillations and therefore presents a background, 
reducing visibility of the oscillations. The contribution 
due to the elastic processes contains both oscillating and 
constant ("elastic background") components, whose rel- 
ative size depends on the relative transmission of the two 
arms of the interferometer. The elastic background is 
minimized when the transmissions through the two arms 
of the devices are of the same order. 

So far we have considered the system of levels with 
strictly alternating parity. Fig. [S^. However, the effect 
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FIG. 6: (color online) Contributions of different levels in AB 
oscillations as functions of bias at zero magnetic flux. 
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FIG. 7: (color online) Coherent (elastic cotunneling) and in- 
coherent (inelastic cotunneling) contributions to the differen- 
tial conductance (shown respectively in solid blue and dashed 
green) . 



of redistribution of the dot population between different 
levels, considered above, can manifest itself in an inter- 
esting way also in cases when the parity of the levels is 
not strictly alternating, such as, e.g., the case shown in 
Fig. [SId, where level 3 has the same parity as level 2. 
The oscillating contribution to the differential conduc- 
tance of such a system is shown in Fig. [51 where the 
calculations were done with the parameters used above, 
except the Coulomb interaction, which was taken to be 
2.4 times bigger. In addition to the first conductance 
switch (which occur when inleastic cotunneling with par- 
ticipation of level 2 begins) there are two more switching 
evens which appear after the second and third inelastic 
onsets (i.e. when level 3 starts being populated from 
level 2 and slows populating of the latter, and when level 
3 starts to be populated directly from level 1). 
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FIG. 8: (color online) Energy diagrams for the two cases con- 
sidered: (a) strictly alternating parity of levels, and (b) levels 
2 and 3 having the same parity. 
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FIG. 9: (color online) The oscillating contribution to the dif- 
ferential conductance (solid blue) and AB phase of conduc- 
tance oscillations (dash green) as functions of bias voltage in 
the case when the parity of levels is not strictly alternating 

(Fig. Eb). 



Thus, the switching mechanism proposed in this paper 
has rather moderate demands in respect to the statistics 
of level parity in a quantum dot. 

One may try to explain the results of Ref. [lOj as origi- 
nating from electrostatic AB effect Indeed, due to their 
small transmission, the quantum dots can be thought of 
as barriers embedded in the two arms of the interferome- 
ter, the whole device thus being an implementation of the 
proposal for measuring electrostatic AB effect made by 
Nazarov-iiii^ The main distinction between the electro- 
static AB effect and the scenario described in this paper 
are the intervals between switchings. In electrostatic AB 
effect the switchings have to occur with a certain bias 
voltage period, determined by the electron time-of-flight 
through the interferometer. This means that switching 
events have to be separated by approximately equal bias 
voltage intervals (the period varying to the extent that 
the density-of-states depends on energy). On the other 
hand, the mechanism proposed in this paper allows for 
the switchings to be separated by unequal bias intervals, 
as indeed is the case in the experimentf^^ Moreover, the 



first switching is clearly correlated with the onset of in- 
elastic cotunneling, which cannot be the case for electro- 
static AB effect. Experimentally these two mechanisms 
can be compared by studying interferometers with the 
same size of quantum dots but different length of the 
arms, i.e. different period of electrostatic AB oscillations. 

Another important way to distinguish between the two 
possible explanations is by observing the characteristic 
behavior of the phase of AB oscillations as it changes 
between and tt. In electrostatic AB effect this change 
is smooth, almost linear in biaSfi^^ii whereas in the ex- 
periment by Sigrist et al. the phase is almost constant, 
except when the bias is close to its switching value, in 
which case the phase flows rapidly between and tt. 
The latter phase behavior can be reproduced within the 
switching mechanism proposed in this paper. However, it 
requires application of many-body scattering theory and 
cannot be done within the one particle scattering ma- 
trix approach employed above, as it incorrectly describes 
phase behavior at small bias. Therefore we leave the dis- 
cussion of these interesting phase phenomena to another 
publicationii^ 



IV. CONCLUSIONS 

We have addressed the findings of the experiment by 
Sigrist et alfiS, We have shown that switching of the phase 
of the AB oscillations of the conductance between and 
TT as a function of bias can be explained as the result of 
transport through the interacting arm of the interferom- 
eter being dominated at different bias by quantum dot 
states of different " parity" . 

The number of switching events is not necessarily equal 
to the number of levels in the interacting arm, which 
agrees with experiment and finds its explanation in com- 
plex redistribution of electron population between differ- 
ent levels. The same redistribution implies that conduc- 
tance switching is not necessarily the sign of alternating 
parity of levels within the quantum dot. 

Finally, the correlation between the first switching of 
the phase and the onset of the inelastic cotunneling, as 
well as unequal separation of the biases at which the 
switching occurs, give reason to think that our expla- 
nation of the experiment is preferable to the one based 
on electrostatic AB effecti ^^i^^ 
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APPENDIX A: SCHRIEFFER- WOLFF TRANSFORMATION FOR A QUANTUM DOT INSERTED IN A 

TIGHT-BINDING CHAIN 

Let us consider the following Hamiltonian, describing a QD sandwiched between two leads, where each level of the 
dot is coupled to a separate channel in the leads. The channels in the leads are assumed to be identical, whereas the 
states in the dot have different energies: 

H = Hl + Hr + Hd + Ht (A1) 



Hd = ^ eafia + ^ X! C*^" = d^da) (A2a) 



2 



-1 



HL=-t «,-lCa,,+C+jC„,j-l) +14 J2 (A2b) 



a,j — — oo Q,j — — oo 

+ 00 +00 



Ht = Y {tLad^Ca.-l + tla^a.-ida + tRad^Ca.! + t*RaCt,lda) (A2d) 

a 

In order to perform Schrieffer- Wolff transformation^'^ we first transform the lead Hamiltonians to the representation, 
in which they are diagonal: 

+0O -1 j-^ 

CRak = y^^<l)k{n)ca,n,CLak = 0fe (n)cQ,„ , 0fe (n) = W - sin (fcn) (A3) 

^ — ' — ' V 71" 

n—1 n— — oo 

The reciprocal relations are 

Ca,« = / dk(j)k{n)cRak, ifn>l;ca,„= / dk4>k{n)cLak, ifn<-I. (A4) 
Jo Jo 

The lead Hamiltonians take form: 

/•TT /"TT 

-ff L = / dk{ek + VL)c'lak'^Lak, HR = y2 dk{ek + VR)c^^^CRak, (A5) 

where 

Efc = -2tcos(A:) (A6) 

is the kinetic energy of electron with quantum number k. 
Now the tunneling Hamiltonian takes form: 

/•TT 

A't = / dk [tLa4'k{~^)d^CLak + t*La4>ki-i)c'^akda + 

a -^0 

tRa4'k{i)daCRak + t*Ra'l>k{i)cR^^.da 1 . (A7) 

In order to perform the SW transformation we need to find an operator 5* that satisfies the equation 

Ht + [S, Hl + Hr + Hd] = 0, (A8) 
then the transformed Hamiltonian will take form: 

H = Hl + Hr + Hd + ^[S, Ht] ■ (A9) 
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We choose operator S as 



F : 



dk 



S = F-F+, 

^RakdaCRak + ^ Biiak,l3'n/3d^CRak + 



-^Lakd^CLak + ^ BLak,l3'ni3d'^CLak 



(AlO) 



This operator is sufficient in the case when the number of channels is two. For a greater number of channels it is 
necessary to add terms that account for the QD being occupied by two, three and more electrons. However, we neglect 

such a possibility and consider only single electron entering or leaving dot. In mathematical terms this is equivalent to 
performing SW transformation an a restricted basis, which includes only the states with zero, one and two electrons 
in the dot. 

Now it is straightforward to show that the coefficients A, B should have forms 



. _ tRa(t)k{l) . _ tLa<t>k{-'i-) 



BLak,f3 = tLa4'k{ — ^) 

The transformed Hamiltonian is 



Ca — Cfe — Vl ■ 
1 



Ea - fife - Vr 
1 

Ca + U -Ck-VR Ca-Ck-VR 
1 1 



(All) 



ea + U -ek-VL ea-ek-Vz 



where 



H = HL + HR + HD + Y,Wj, 



Wl = ly2 [ dk [ARc,kt*RM^) + ALaktlM-'^) + 

A*RaktRa(Pk{l) + Al^f,tLa(l)ki-'i-))na = ^ Ae^na 



(A12) 



(A13) 



describes the correction to the dot energies due to coupling to the leads. In the following we will neglect this term 
the appropriate corrections can be included in energies €(x of the Hciniiltonicin Ho- 



ct 

1 



(A14) 



itRat'ha 



~^ Ca-ep - Vr 
1 1 



ep - Vr 



cf-Vl 



Ca-ep- Vl 



W2 describes clastic cotunncling through the QD without account for interactions. Following SW we replaced in the 
denominators electron kinetic energy, Ck, by the Fermi energy. It is essential to keep in mind that, since the leads are 
kept at different bias, their chemical potentials are different: Hl = + Vl, Mij = cf + Vr. 

W3 = - V / dk dk' BRak,f3tRi3(l)k'{i)d0d+CRakCR0k> + 
^ „, n-/-^ JO Jo L 



BRak,l3tL04>k'{-'^)d'pd'^CRakCLI3k' + BLak,l3tRl34>k' {^)d'^ d^CLakCRfjk' + 

BLak,l3tL04>k'{^)d'^d^CLakCLl3k' + H. 



(A15) 
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describes simultaneous tunneling of two electrons into or out of the dot. In the following we neglect this term, 
since the corresponding processes are energetically unavailable. 



W4= ^ \^v'^r<^'li^ld'^di3Ca,l+vf"c'^_^d'^di3Ca,l + 



vtc+_,dtd0C^,.A , (A16) 



'^rr — -^^RatRf:} ( " TTT 77 ~ 1 77 — ^ 



Bex I 1^ I - - I 
% = l^tLatL/S [ — -r — + 



2 



ea + U 


- ep - 


Vr Coc - 


ep 


-Vr 




1 




1 




e0 + U 


-ep- 


Vr ep - 


ep 


-Vr 




1 




1 




ea + U 


- ep - 


Vl Ca - 


ep 


-Vl 




1 




1 




ea + U 


- (F - 


Vl ep- 


fF 


-Vl 




1 




1 




ea + U 


-€f- 


Vl Ea - 


ep 


-Vl 




1 




1 




ep + U 


- ep - 


Vr ep - 


ep 


-Vr 



^LatRp I - — FT \ T7 \ \ ~ ^ 



W4 describes inelastic cotunncling through the quantum dot, i.e. the tunneling events when the dot changes its state, 
while the incident electron is transferred from one lead to the other or reflected back. 



W5 = - / dk {BRak,f3t}ia(l'k{^) + BLak,/3t*La(l>k{-^) + C.C.) UaUp = 

^ AUcprianp (A17) 

W5 describes correction to the Coulomb energy and is expected to be small; in the following we neglect this term. 
Finally, we also have 

We = - J2 (C^r^/JC+iCa,! + <"n/3< _ic„,i + <;«n0C+ ic„,_i + vrnpcl_,ca,-i) (A18) 

1 1 



vr = \tLc 



2 



ea + U -ep - Vr ea- ep - Vr) ' 
1 1 



ea + U - ep - Vl ea - ep - Vl 



- ^tRatl^ {^^ + u-ep- Vr ea-ep-VR^ 
1 1 



ea + U - ep - Vl ea - ep - Vl 

a< 



aa _ („,aa\* 



which describes the correction to elastic cotunncling contribution, W2, which arises when the dot is occupied by an 
electron {W2 is clastic cotunncling through an empty dot: "hole tunneling"). 

For the reference arm, in which we neglect the Coulomb interaction, the Schrieffer- Wolff transformation readily 
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produces the only tunneling term: 

^'■e/ = ~ X! i'^rrC^^iCa,! + UJlrC^-lCa,! + UJrlC^.iCa-l + LOuC^ _-^Ca-l) (Al9) 
a 



eo - CF - Vr £0 - - Vl 

where eo and tj^ Q^tn^Q are the level energy and the tunneling amplitudes for this arm, which we assume to be 
independent on the channel index. 

We neglect Wi_5 describing corrections to the coefScients and W3 describing two-particle tunneling, but keep the 
terms W2,4.6 which, together with H^, H^, Ho and Wref form Hamiltonian Ti., Eq. ((T|), used in this paper. 



APPENDIX B: SCATTERING MATRIX 



Using the scattering wave function in the form of Eq. (jj]) , it is easy to show that the amplitudes satisfy equations 



M 

M = 











Mll 


Mlr 


Mrl 


Mrr 



Br 



,^ S = -M-^M: 



Mll Mlr 
Mrl Mrr 



The only non-zero elements of matrices , M are 

{vZ. + {l-6^^p)vT)e''',MLR 



Mfr'"^ = t - [vfi + (1 - 5^^p)vfne'''\M^^-'^'' 



LL 



'■LL 

MfR"" 

Mr^'^' = - + (1 - 5«.^)<^)e^^^^lSi^''" 



= (1 - &a,p)Vii 
= (1 - <5a,/3)W/"' 
= (1 - ^a,/3X 



a!3 ^ika 



and 



Mlr""" 



t - Ki + 



Uf^-P = t - Kr + (1 - 6^,f,Kn e-'''^,Mf/'^ 

I 



(1- 
(1- 
(1- 
^(1 



a/3 



(Bl) 
(B2) 



(B3) 



(B4) 



APPENDIX C: PARAMETERS USED IN THE 
RESULTS SECTION 

The level energies were taken to be ei = —AmeV, ti = 
— .ZhmeV , tz = —■25meV; the tunneling amplitudes are: 
tL,i = 50iieV,tL.2 = 88.9fieV,tL,3 = 150/iey, tL,i = 
50iJ,eV,tL,2 = — 88.9 fj.eV,tL ^3 — IbO^ieV, which corre- 



spond to the level widths: T^ i — Tr^i ~ .25fieV, T^ ^ = 
Tr^2 = ■79iJ,eV,TL^3 = Tr^s = 2.25^6^ (in tight-bonding 
picture level width is defined as F = |i^,aP/i)- The 
Coulomb interaction is U — bmeV. The parameters 
of the reference arm are: eg = —ImeV, iL,o = tR,a = 
200^eV, Tl.o — Fl.o = AfieV. The temperature is taken 
to be T = 60mK. 
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